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1. Introduction 

Recent experiments on the spin Hall effect PU 0] and the spin galvanic and circular 
photogalvanic effects 0] have stimulated renewed interest in the behavior of electrons 
and holes in semiconductor nanostructures in an external magnetic field. Such 
experiments are often modeled using an effective-mass or envelope-function Hamiltonian 
E] that may include a generalized Rashba coupling jTj to account for spin-splitting 
effects. The magnetic field is included via the minimal coupling substitution k — > 
k + A/c, which is well justified for bulk semiconductors [HJ I§1 HUj. 

However, Mlinar et al have recently demonstrated that the results of 

such calculations can depend strongly on the particular choice of envelope-function 
Hamiltonian. The standard "symmetrized" Hamiltonian J2| applied to GaAs/AlAs 
and InAs/GaAs quantum dots in a strong magnetic field yielded energy eigenstates 
with arguably unphysical properties [JJ] — in contrast with the physically reasonable 
predictions ^Tj obtained from an asymmetric Hamiltonian fHJ Ej • This may be viewed 
as a modern confirmation of Luttinger's original work jU] demonstrating the importance 
of including antisymmetric terms in the Hamiltonian. 

Nevertheless, it is not clear that the model used by Mlinar et al [II] provides a 
correct description of the heterojunction regions. The antisymmetric part of the valence- 
band Hamiltonian derived in reference [13] is in fact just the valence-band Rashba 
coupling [7] that is of interest for spintronics applications |T5] . However, a recent 
first-principles envelope-function theory for lattice-matched heterostructures ^B] has 
established the existence of several other interface terms of the same order of magnitude 
as the Rashba term. These terms are not included in standard effective-mass models 
for nanostructures, although some authors |I7J have included phenomenological surface 
terms in their study of the magnetic field problem. 

The purpose of this paper is to establish a rigorous theoretical foundation for the 
effective-mass theory of semiconductor nanostructures in a magnetic field. This is done 
by extending the first-principles multi-band effective-mass Hamiltonian of reference ^E] 
to include the effects of a slowly varying static external magnetic field. In this theory 
[To] , the nanostructure pseudopotential is treated as a perturbation of a bulk reference 
crystal (e.g., Alo.5Gao.5As for a GaAs/AlAs quantum dot), with the self-consistent linear 
and quadratic response terms included in the framework of the k ■ p perturbation theory 
of Leibler [IB]. 

Until recently, such a first-principles analysis would not have been possible, due 
to the lack of a generally accepted Hamiltonian to describe the coupling between the 
magnetic field and the nonlocal part of the potential energy. Such nonlocal terms arise 
for two reasons. At the fundamental level, the self-energy operator (which accounts for 
electron-electron Coulomb interactions) in Dyson's equation is nonlocal. At the practical 
level, the pseudopotentials used in most first-principles calculations (which account for 
the influence of core electrons) are also nonlocal. 

Pickard and Mauri ^H] have solved the second problem by deriving a coupling 
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for norm-conserving pseudopotentials in the local density approximation (LDA) [20] 
that reproduces the results of all-electron LDA calculations. Ismail-Beigi et al [21] 
have proposed a more general coupling for arbitrary nonlocal potentials, but since this 
coupling, in the special case of norm-conserving pseudopotentials, does not agree with 
all-electron calculations ^H], its validity in the general case is questionable. 

The present work is therefore limited to the case of norm-conserving 
pseudopotentials within LDA. Although LDA calculations do not accurately predict 
energy gaps in semiconductors, this at least provides a working model that is self- 
consistent and accounts for the true atomic structure of a heterointerface. The results 
obtained here should be representative of the general features obtained in a more 
complete theory that includes the nonlocal quasiparticle self-energy. 

The results derived in this paper confirm the existence of interface terms in the 
Hamiltonian that were not included in the calculations of Mlinar et al ^T] . An explicit 
numerical calculation of these terms is in development and will be reported elsewhere. 
However, the present analytical results show that these terms are of the same order of 
magnitude as the valence-band Rashba terms studied in reference [TT], which suggests 
that they may have a comparable influence on the energy eigenvalues and eigenstates 
in an external magnetic field. 

In addition, an analytical demonstration of the difference between the coupling 
Hamiltonians of Pickard and Mauri [19 and Ismail-Beigi et al [21] is established for 
the first time. In particular, the nonlocal pseudopotential contributes to the effective 
magnetic dipole moment of the electron in the first case, while in the second case it does 
not. 

The paper begins in section [2] with the derivation of exact expressions 
for the magnetic-field coupling Hamiltonian in coordinate and momentum space. 
Approximations suitable for slowly varying envelope functions are developed in section 
E] based on a power series expansion of the nonlocal pseudopotential in momentum 
space. These results are used to construct an envelope-function Hamiltonian in section 
E] The contributions derived from k ■ p perturbation theory are added in section 
where the complete expression for the effective-mass Hamiltonian is given. The results 
are discussed in the concluding section El Throughout the paper, Hartree atomic units 
(h = m = e = 1) are used. 

2. General form of the coupling Hamiltonian 

The magnetic coupling Hamiltonian derived by Pickard and Mauri JH] is applicable 
to nonlocal potentials formed by a superposition of localized nonoverlapping norm- 
conserving pseudopotentials; i.e., 



R 

where Vf^(x, x') is the nonlocal part of the pseudopotential for an ion at position 
R. In this case, the potential in the presence of a slowly varying magnetic field 




(1) 
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B{x) = V x A(x) is given by [TH] 

V*{x, x') = V Vg{x, x') exp ( - [ A{x") ■ dx") , (2) 

R \ c Jx^R-^x' J 

where x — > R denotes an integral over the straight-line path from x to R. If A(x) is 
chosen to be periodic [TTH ITTH 12*2*] , it can be expanded in a Fourier series, which permits 
the line integral to be evaluated as 

/ A ■ dx" = (x'-R)-Y] A(k)e ik - R f[ik ■ (x' - R)], (3) 



in which 



2 

X X 



/W = _ = i + _ + _ + .... (4) 

Since the coupling (0) is valid only for slowly varying B(x) it is assumed that A(k) 
is significant only for small k (i.e., where a is the lattice constant). The power 

series in (0J) therefore converges rapidly, because the range of the nonlocal potential (i.e., 
the core radius) is less than a. From this power series one obtains 

f A- dx" = (x> - R) ■ A(R) + ±(a£ - R x )(^ - 

+ - Rx)K - R,){< - R »)^^ u + ■ ■ ■ > (5) 

in which a sum over repeated Cartesian coordinate indices A, /i, v is assumed, and 
dAx/dRp = (dA x /dx^)\ x=R . Therefore 



A ■ dx" = (x' - x) ■ A(R) 

x^R-^x' 

+ ^[{A - Rx)K - *m) - (** ~ R ^ - ^)lf^ + ■■■■ ( 6 ) 

Note that equations (jSJ) and (jSJ depend only on the symmetric part of the matrix 
dAx/dRp, i.e., they contain no terms proportional to B or its derivatives. As shown in 
the Appendix, this is generally not true for other coupling Hamiltonians. 

For the potential energy of a lattice-matched heterostructure, it is convenient to 
use a more explicit notation: 

V(x,x') = ^2f%v a (x- R a ,x' - R a ), (7) 

a,R 

in which R is now a Bravais lattice vector (for the periodic bulk reference crystal), a 
labels an atom at position r a in the unit cell, and R a = R + T a . The function /j| is 
the atomic mole fraction of atom a at position R a , while v a (x, x 1 ) is the nonlocal ionic 
pseudopotential for atom a. 

With the coupling to the magnetic field given by Eqs. (0) and (|SjL the Fourier 
transform of the field- dependent potential can be written as 

V(k, k') = f^ k - k '> Ra exp^Or, x')/c]v a (K, K'), (8) 

a,R 
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in which K is the kinetic momentum 

K = k + ^-A(R a ), K' = k' + ^A(R a ). (9) 
The operator function £ is defined by 

£(x, a; j = ^(a^a^ ~~ x x x n) Qj^ a 3! ^ >?n Xv ~ x ^ x ^ x v) Qj^aQj^a ~<~ ' ' ' > \^) 

in which x\ = id/dk\ and x'^ = —id/dk'^. These expressions are used to develop an 
effective-mass Hamiltonian in what follows. 



3. Approximations to the Hamiltonian 

Effective-mass theory is based on approximations valid for k cnv a <C 1, where k env is a 
typical wave number for a slowly varying envelope function. This condition will hold, for 
example, in a wide quantum well of width L ^> a, for which fc env ~ 2n/L. The external 
magnetic field B must also be weak enough for the magnetic length L B = yc[B to 
satisfy L B ^> a [TU] . 

In addition to these standard conditions, if K is to be treated as a small quantity of 
the same order as A;, the vector potential must be small enough to satisfy aA/c <C 1. For 
a sinusoidal field with B ~ k B A [22], where k B is the wave number of the external field, 
this condition is equivalent to k B L 2 B ^> a. Combining this with the previous condition 
L B ^> a, the new condition will be satisfied if k B L B > 1. Thus, very small values of 
k B place an additional constraint on L B . This constraint is not usually mentioned jlUj . 
but it is implicit in any effective-mass treatment where K is assumed to be of the same 
order as k. Therefore, effective-mass calculations for truly constant B fields (which are 
obtained in the limit k B — > 22J can only yield results that are valid asymptotically 
in the limit B — > 0. The asymptotic nature of effective-mass theory for constant B was 
suggested by Kohn [23], and stated more explicitly by Foldy and Wouthuysen [21] in 
connection with the nonrelativistic approximation to the Dirac equation. 

In the Luttinger-Kohn representation \nk) [H], the wave function (nk\ip) is just a 
unitary transformation of the momentum-space wave function (k + G\ip) |25| : 

(nk\^ = Y,Ul G (k + G\^). (11) 

G 

Here G is a reciprocal lattice vector of the reference crystal, and the spinor U n a 
is a Fourier series coefficient for the n th zone-center Bloch function of the reference 
crystal. For slowly varying envelope functions F n (x), the Fourier transform F n (k) = 
(nk\ip) (and hence {k + G\ip)) is significant only for small k. (A review of the 
numerical evidence supporting the use of the slowly varying envelope approximation 
in abrupt heterostructures is given in reference [16|.) Thus, slowly varying envelope 
functions probe the nonlocal pseudopotential only within a small neighborhood of the 
reciprocal lattice vectors G. Within this neighborhood, the ionic pseudopotential can 
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be approximated by the leading terms in its Taylor series expansion [To] : 



v a (k + G,k' + G>) = w g£ + k x v^, a + k' x v^, a + k x k^ a + k' x k'^, 

+ k x k'/*S a + ..., (12) 

where Vqq, is an expansion coefficient. This expansion is valid because v a (x,x') is 
localized, hence v a (k, k') is analytic. For simplicity, terms beyond the second order are 
not written out explicitly here; however, the present effective-mass theory [TH] retains 
terms of order k 4 in the Hamiltonian of the reference crystal [see equation ()17j) below]. 

By definition, the expansion coefficients are symmetric with respect to arbitrary 
permutations of the indices on one side of the vertical line — for example, Vqq, = 
v GG^ a ■ However, Vqq, ^ v Gc> a - F° r indices on opposite sides of the vertical line, 
the only symmetries are those obtained from hermiticity and time reversal, such as 

With the above restrictions on the magnetic field, K is a small quantity of the 
same order as k, and the potential (jHJ) can be expanded in a Taylor series about the 
reciprocal lattice vectors G. This yields 



,(A|> 



(•|A)a 



,(A/4)a 



(■|A/i)a 



V(k + G,k' + G') = J2f^ 



-i(k-k')-R a -i(G-G')-T° 



a,R 



U GG' + A A^GG' + JX X U GG' 



(A[-)a 



,(-|A)a 



+ ( K X K, + 



id Ay 



A \„,(Am|-)« 



cdK* 



v 



GG' 



■ (/ v ; a;, 1^Hv a 



GG' 



> L 'GG 



+ 



(13) 



where the expansion coefficients are the same as in ([12)1 . and the terms involving 
dA x /dR^ are the leading contributions from (jl(Jj) . 

4. Envelope-function Hamiltonian 



These results are now used to construct an envelope-function Hamiltonian valid for 
slowly varying envelopes. In the Luttinger-Kohn representation, the matrix elements of 
V are given by the unitary transformation 

V nn ,{k, k') = U ] nG V{k + G,k' + G')EW (14) 



G,G 



When the expansion Ijlflj) is inserted into equation (|14jl . one obtains transformed 
expansion coefficients such as 

rt jXn\-)a TT ._.„-i(G-G')-r a 



y{\n\-)a 
nn' 



G.G' 



in which iV is the number of unit cells in the crystal. It is also convenient to introduce 
coordinate-space matrix elements of the form 

^ |,) (*) = E^ l ' )0 /°(*), (16) 
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in which f a (x) is a macroscopic average (see reference (THj) of the discrete atomic 
distribution function /«. Note that Vffl'\x) = V^'\x) = [V^\x)\* and 
V^"\x) = [V^\x)Y . 

With these notational conventions, the Fourier transform of the matrix element 
(|T3J) yields the following contribution from the nonlocal pseudopotential to the envelope- 
function Hamiltonian: 

V n Ax,p) = V$\x) + PxV^\x) + V^\x)P x + P x P,V^\x) + V^\x)P x P, 

+ P x V^\x)P a + V r x jP x P a P u + V n x n rP x P,P K P u , (17) 

in which the canonical and kinetic momentum operators are p = — iV and P = 
p + A(x)/c. The derivation of this result makes use of the usual approximations JU] 
for slowly varying potentials A(x) and envelope functions F n (x). In other words, it was 
assumed that since only small values of k are of interest, the error involved in replacing 
the Fourier transform of the discrete product f^A(R a ) with the Fourier transform of the 
continuous product f a (x)A(x) is negligible, and the upper limit of any convolutions in 
k space can be extended to infinity. This may be called the local approximation because 
it reduces equation (jl7j) to the form of a local differential operator; see reference ^Hj for 
a discussion of the error involved in this approximation. 

In keeping with the perturbation theory of reference jTHj, the full position 
dependence of f^ was retained in the terms quadratic in P in equation p7jl . However, 
in the cubic and quartic terms, /j^ was replaced by the atomic distribution function f a of 
the reference crystal (i.e., the heterostructure perturbation 9% = f^—f 01 was neglected). 
Thus, the coefficients of the cubic and quartic terms are given by the constants 

yXjJ,v _ y(Xflv\-) , y(Xl*\v) _|_ yWfU>) _|_ y{-\Xia>) /io\ 
nn' nn' ' nn 1 ' nri' ' nn 1 ' \ / 

yXflKU _ y(\HK,v\-) y{\^K\u) y(\fJ,\Kv) y(X\flKU) y(-\X^KU) , qX 

Tin' nn 1 ' nn' ' nn' ' nn' ' nn' ' V / 

in which the various terms are just expansion coefficients that were omitted in equation 

Because of the general lack of symmetry of the expansion coefficients in equation 
(|12)1. the noncommuting operators in equation (j!7)) have a definite operator ordering. 
The term involving V^,^ (x) has the same operator ordering as the Hamiltonian studied 
in the work of Mlinar et al [TT] . and includes terms antisymmetric in A and /i. However, 
that Hamiltonian |T51 ITi] does not include any terms with the operator ordering of 
V^\x) and V^\x). 

As might be expected from the work of Luttinger equation (fTTj) is not the same 
as what would be obtained by substituting p — > P into the Hamiltonian of reference 
[TT)j . The position-dependent terms describing the linear response to the heterostructure 
perturbation do have this form, since care was taken in reference jTB] to retain the 
asymmetric parts of these material parameters. However, the bulk coefficients derived 
here [such as (fTKj). (JTOJ) , and the bulk part of V„^J^] are asymmetric in the Cartesian 
indices, whereas those obtained from the expansion of v^, G ,(k) = v a (k + G,k + G') 
in reference JIE] were symmetric. These asymmetric coefficients generate terms in the 
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Hamiltonian proportional to the magnetic field, such as the dipole term in equation (|2(Jj) 
below. 

5. Other contributions to the envelope-function Hamiltonian 

Within LDA, the self-energy of the electron is a local potential that does not couple to a 
static magnetic field [T^]. Thus, the only other term in the Hamiltonian that depends on 
the magnetic field is the kinetic energy \p 2 ■ To account for the intrinsic dipole moment 
of the electron, this can be written as |(<r • p) 2 , where er is the Pauli matrix vector; the 
minimal substitution p — > P then generates the electron g factor go = 2 automatically 

In the Luttinger-Kohn representation, the kinetic energy produces k ■ p terms of the 
usual form [HI EU E3 a finite-dimensional multi-band effective-mass equation can then 
be derived using the standard canonical transformation method of perturbation theory 
[3 El EEH UE] , with the k ■ p interaction and the heterostructure potential treated as 
perturbations [TT)| I18j. The details are essentially the same as those given in reference 
|16j . which will not be repeated here. The outcome is a Hamiltonian with the same 
qualitative structure as that shown in equation (fTTj) . but with all of the coefficients 
renormalized [e.g., V^^'\x) — > V^^'\x)\ to account for perturbative corrections. 

This Hamiltonian can then be rearranged into a more convenient form |18j . which 
yields the effective-mass Hamiltonian 

H nn ,(x,p) = EJ nn , + {P x , ^ n ,(x)} + {{PxP,}, D^\x)} - ti n ,(x)B x (x) 

+ P X P a P K C X S + PxPaP.PuQnT + W nn ,{x) + V ( X )5 nn , 

+ d x Z x n ,{x) + d x d u Y n x :,(x) + {P Xl d^ n ,(x)} + [d^,(x)]P x 

+ Z^dMx) + Y^,d x d^{x) + t%,{P x , dM*)} + Kn\ d M*)\P^) 

in which {A, B} = {AB} = ^(AB + BA) is the symmetrized product, and d x = d/dx x 
acts only on the function immediately to its right. The various coefficients in this 
equation are the same as those defined in reference [15] . with some modifications to be 
discussed below. Therefore, the definitions given in ^5] are not repeated here. 

The first term E n is just the zone-center energy of the reference crystal, while 
7r^ n ,(x) is the position-dependent momentum matrix of the heterostructure. The 
Luttinger D matrix [2j, which is half the inverse effective mass tensor, is defined by 

D X M = l(Vnn> + iej^,) + V%f-\x) + V^\x) + V^\x), (21) 

in which e Xnu is the antisymmetric unit tensor. Functions such as V^,(x) have the 
same form as V^f (x), but include additional k ■ p renormalization terms that are 
given in ^5]- However, v}t^'\ x ) * s no l° n ger symmetric: V^f (x) ( x )- 

The symmetric part of the tensor (|21j) is written as D^, (ac), while the 
antisymmetric part determines the effective magnetic dipole moment 

/&'(*) = h^ D nn>i X )- (22) 
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The symmetric terms V^i (x) and V^,(x) do not contribute to (}22|) : the sole 
contribution from the nonlocal potential is the asymmetric term V^j (x). Note, 
however, that this result depends on the initial choice of (0) as the magnetic coupling 
Hamiltonian. As shown in the Appendix, for the generalized Peierls coupling of Ismail- 
Beigi et al [21], the nonlocal potential contributes nothing to the dipole moment (|2*2*Jl . 

The cubic (C n ^f) and quartic (Q n ^ U ) dispersion terms are just renormalized 
versions of (fTHj) and (fT^j) . The term W nn >(x) represents linear and quadratic 
contributions to the effective potential energy, one portion of which is V^,\x) from 
equation (fTTj) . The screened potential (p(x) accounts for long-range multipole potentials 
generated by the electron-electron Coulomb interaction [T75] . 

The terms Z* n ,(x), Y n £,(x), T n ^,(x), and $^,(a;) are short-range interface terms, 
the contribution from which vanishes in bulk material because the spatial derivatives 
d\ are zero there. Z generates a 5-like interface mixing of light and heavy holes |27j . 
while Y is similar to the Darwin term in the nonrelativistic approximation to the Dirac 
equation. The term $, which is antisymmetric, is a generalized Rashba coupling 7\ that 
can be viewed as an effective spin-orbit interaction. This is precisely the antisymmetric 
term that was examined in the work of Mlinar et al jTIJ. The T term has a similar 
structure (i.e., an interface term that is linear in P) but is symmetric. 

The contributions from Z, Y, f , and $ are similar to those from Z, Y, T, and 
but since they involve the long-range multipole potential f(x), they are not as well 
localized at the interface ^B]. These terms are not included in most envelope-function 
theories. However, in a quantum well or other two-dimensional system, these terms have 
no qualitative effect, as they merely renormalize the short-range contributions fB]. This 
is not strictly true in quantum wires or dots, but even in these cases it may be possible to 
replace the long-range contributions with renormalized short-range terms, since Zunger 
et al |2~%| EHJ ED) have shown that LDA wave functions in a variety of nanostructures 
can be accurately reproduced using carefully fitted empirical pseudopotentials, which 
by definition have no long-range multipole terms. 

6. Conclusions and discussion 

In this paper, a first-principles multi-band effective-mass Hamiltonian was derived for 
lattice-matched semiconductor nanostructures in a slowly varying magnetic field. The 
theory applies to systems described by norm- conserving nonlocal pseudopotentials in 
the local density approximation, and is based on the magnetic coupling Hamiltonian 
derived by Pickard and Mauri [TH]- The main result, given in equation (|2()jl . is similar 
to what would be obtained by applying the minimal substitution p — > p + A/c to 
the Hamiltonian derived in reference |16j . However, there were also several additional 
terms, derived from asymmetric matrix elements of the nonlocal pseudopotential, that 
were not included in reference [HI]- As shown in the Appendix, these terms are controlled 
by the choice of magnetic coupling for the nonlocal pseudopotential. The present work 
provides the first direct analytical demonstration of the difference between the magnetic 
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coupling Hamiltonians of Pickard and Mauri [F3\ and Ismail-Beigi et al [21] . Namely, in 
the former case the nonlocal potential contributes to the effective magnetic dipole 
moment, while in the latter case [21] it does not. 

Turning now to the significance of the various interface terms in the Hamiltonian 
(ffilj) . the work of Mlinar et al [TT] has established the importance of including the 
generalized Rashba coupling ^n'( x )- m their calculations, the magnitude of this term 
was estimated using the expressions given in references J3] and [Ej, which are equivalent 
to the assumption that $^,(cc) is the same as the antisymmetric part iD£$(x) of the 
effective mass tensor (jUJ), or that $^,(cc) = ce arv \i v nn i{x) . (Note that for T 8 states in 
zinc-blende, /Ji u (x) = [k(x)J u + q(x)J^]/c, where J is an angular momentum | matrix, 
and k and q are Luttinger parameters [§]. The calculations in reference ^JJ included 
only the contribution from k.) As discussed in reference ^Hj, this assumption provides 
a reasonable order-of-magnitude estimate, but cannot be relied upon for quantitative 
accuracy. Thus, the large qualitative impact of $^,(ai) discovered by Mlinar et al [TT] 
is almost certainly a real physical effect, even though the quantitative details of their 
calculation may not be completely accurate. 

In addition to the Hamiltonian pUj) also includes the interface terms 

Z^ n ,(x), Y^,(x), and r^,(:r), which were not studied by Mlinar et al [TT] . A 
quantitative calculation of these matrix elements is now in progress (and will be reported 
elsewhere), but a rough idea of their significance can be obtained from dimensional and 
symmetry analysis. The terms Z* n ,(x) and Y^,(x) affect the boundary conditions on the 
envelope functions, but do not contain any direct contribution from the magnetic field. 
Therefore, their influence on the magnetic-field-dependent properties of nanostructures 
is probably not as strong as that of $^,(:r). Nevertheless, their contribution may still 
be significant — especially the ^-function term Z^ n ,(x) [21], which is of lower order [TE] 
than the other interface terms. 

On the other hand, the term involving T^^x), like $ n ^, (x), is linear in the vector 
potential. The symmetry analysis of reference JH] shows that this term is not present in 
the Tq conduction band of zinc-blende materials, and that in the Tg valence band it is a 
relativistic effect that vanishes when spin-orbit coupling is neglected. Its magnitude will 
therefore be smaller than that of $^,(:r), but it may be important in the description of 
spin-splitting effects. 

Therefore, even though the calculations of Mlinar et al jTJJ have included what is 
probably the dominant magnetic-field-dependent interface term, the remaining terms 
are not likely to be negligible. The detailed study of these terms, however, is left to 
future numerical work. 
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Appendix 

This appendix examines the consequences of replacing the coupling Hamiltonian of 
Pickard and Mauri [TH] with that of Ismail-Beigi et al |2T| . who have proposed a 
generalization of the well known Peierls phase [HU E2j : 



V n \x, x') = V n \x, x') exp ^ J A{x") ■ dx"^j . (A.l) 

Unlike equation (j2J), this does not rely on separating V nl (x,x') into a sum of localized 
atomic potentials — although, for the application considered here, it is assumed that the 
potential has this form. The line integral can be evaluated by replacing R with x in 
equation (jSJ) and then expanding A(x) in a Taylor series about x = R, which gives 

/ A ■ dx" = (x' - x) ■ A(R) + h{x' x - R x )(x> - R p ) - {x x - R x )(x, - R P )}^ 

Jx^x' 1 aK p 

-i(, A - BA )«- B j(^-^) + .... (A.2) 

The last term, which was not present in equation (JEJ), is proportional to the magnetic 
field. 

If the derivation in sections E2 and 0] is now carried through in the same way as 
before, one finds that the last term shown in equation ()13|) is modified as follows: 

KxK'/Sg" - ( K ^K + ^B v {R a )) v ( G l $ a . (A.3) 

The effect of the additional term is to cancel the contribution from V^^\x) to the 
effective dipole moment ((221) • Thus, for this coupling, the nonlocal potential makes no 
contribution to the dipole moment. A similar conclusion was reached by Kane [HH]; who 
used the Peierls coupling (jA.lj) for the special case of a constant magnetic field. 

Therefore, different choices of coupling between the nonlocal potential and the 
magnetic field have a direct impact on the effective dipole moment of the electron. 
This provides an explicit analytical demonstration of the differences demonstrated 
numerically in reference ITTJj . 
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